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Abstract: There is a growing amount of evidence that QCD (and four-dimensional 
gauge theories in general) possess a hidden symmetry which does not exhibit itself as 
a symmetry of classical Lagrangians but is only revealed on the quantum level. In this 
review we consider the scale dependence of local gauge invariant operators and high- 
energy (Regge) behavior of scattering amplitudes to explain that the effective QCD 
dynamics in both cases is described by completely integrable systems that prove to be 
related to the celebrated Heisenberg spin chain and its generalizations. 
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1 Introduction 



QCD is a four-dimensional gauge theory describing strong interaction of quarks and 
gluons. There is a growing amount of evidence that QCD (and Yang-Mills theories in 
general) possess a hidden symmetry. This symmetry has a dynamical origin in the sense 
that it is not seen at the level of classical Lagrangian and manifests itself at quantum 
level through remarkable integrability properties of effective dynamics. 

The simplest example which allows us to explain integrability phenomenon is a 
process of deeply inelastic scattering (DIS) of an energetic hadron off virtual photon, 
7*(g) + h{p) — > everything. This process played a distinguished role in early days of 
QCD development and it led, in particular, to important discoveries such as QCD fac- 
torization and formulation of parton model for hard processes (see e.g. |1|). The total 
cross-section of of DIS process is related by the optical theorem to imaginary part of the 
forward scatteting amplitude j*(q) + h{p) — > 7*(g) + h{p) (see Fig. 1). It is parame- 
terized by the so-called structure functions F(x, q 2 ) depending on the photon virtuality 
q 2 < and dimensionless Bjorken variable < x < 1. The latter is related to the total 
center-of-mass energy of the process as s = {p + q) 2 = —q 2 (l — x)/x. 



Figure 1: The total cross-section of deep inelastic scattering j*(q) + h(p) — >• 
everything is related by the optical theorem to imaginary part of the forward 
scatteting amplitude. Solid and wavy lines denote quarks and gluons, respec- 
tively. 

The integrability has been first discovered in Refs. |2-4| in the study of high-energy, 
s ^> — q 2 (or equivalently x — > 0) asymptotics of F(x,q 2 ). Experimental data in- 
dicate that the structure functions increase in this limit as a power of the energy, 
F(x,q 2 ) ~ (l/a;) w , in a quantitative agreement with the Regge theory prediction. At 
weak coupling, the same behavior can be obtained through resummation of perturba- 
tive corrections to the structure functions enhanced by logarithm of the energy |5|. The 
structure functions obtained in this way satisfy nontrivial multi-particle Bethe-Salpeter 
like evolution equations |6,7|. These equations have resisted analytical solution but a 
breakthrough occurred after it was found |2-4 that, in multi-color limit, these equations 
can be mapped into a Schodinger equation for a completely integrable quantum (non- 
compact) Heisenberg SX(2,C) spin chain. This opened up the possibility of applying 
the quantum inverse scattering methods for the construction of the exact solution to the 
evolution equation in planar QCD. 

Later, similar integrable structures have been found in Refs. |8-10| in the study of 
dependence of the structure functions F(x, q 2 ) on the momentum transfered q 2 . At large 
Q 2 = —q 2 , the operator product expansion can be applied to expand the moments of the 
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structure functions in powers of a hard scale 1/Q 



fdxx^F^a 2 ) = CNAa n f 2)) ( P \0 N , L \p)^ ■ (1.1) 
Jo L > 2 V 

Here the expansion runs over local composite gauge invariant operators ( Wilson opera- 
tors) of Lorentz spin N and twist L. The corresponding coefficient functions cjv,i(a s (Q 2 )) 
can be computed at weak coupling as a series in the QCD coupling constant a s (n 2 ) = 
g 2 /(Air) normalized at /i 2 = Q 2 . At the same time, the matrix element of the Wilson 
operator with respect to hadron state (p\On,l\p) h 2 =q 2 is a nonperturbative quantity. Its 
absolute value can not be computed perturbatively whereas its dependence on the hard 
scale Q 2 is governed by the renormalization group (Callan-Symanzik) equations 

/^(plOj&ip) = -1 { nU*sMO { «1\p) . (1.2) 

Here we introduced the superscript (a) to indicate that for given N and L there are 
a few Wilson operators parameterized by the index a. The Callan-Symanzik equation 
(1.2) has the meaning of a conformal Ward identity for the Wilson operators with the 
anomalous dimension 7^*^(a s ) being the eigenvalue of the QCD dilatation operator (see 
e.g. review |11|). 

The Wilson operators are built in QCD from elementary quark and gluon fields and 
from an arbitrary number of covariant derivatives. In general, such operators mix under 
renormalization with other operators carrying the same Lorentz spin and twist. Diag- 
onalizing the corresponding mixing matrix we can find the spectrum of the anomalous 
dimensions 7^^(a s ). For the Wilson operators of the lowest twist, L = 2, the anomalous 
dimensions can be obtained in the closed form |12|, whereas for higher twist operators 
the problem becomes extremely nontrivial already at one loop due to a complicated form 
of the mixing matrix |13|. Quite remarkably, the spectrum of the anomalous dimensions 
can be found exactly in QCD in the sector of the so-called maximal-helicity Wilson op- 
erators. The reason for this is that the one-loop mixing matrix in QCD in this sector can 
be mapped in the multi-color limit into a Hamiltonian of the Heisenberg SL(2,M) spin 
chain |8-10|. The twist of the Wilson operator L determines the length of the spin chain 
while the spin operators in the each site are defined by the generators of the 'collinear' 
SX(2,R) subgroup of the full conformal group |14, 15|. As a result, the exact spectrum 
of one- loop anomalous dimensions can be computed with a help of Bethe Ansatz |16|. 

Let us now examine the relation (1.1) for large Lorentz spin, N ^> 1. This limit 
has important phenomenological applications in QCD |17,18|. It is known |12| that the 
anomalous dimensions of Wilson operators grow as their Lorentz spin increases. As a 
consequence, the dominant contribution to (1.1) only comes from the operators with 
the minimal anomalous dimension 7^ L = mine 7^. Quite remarkably, this anomalous 
dimension has a universal (twist L independent) logarithmic scaling behavior at large N 
to all loops 1 19, 20] 

^ L = 2T cusp (a s ) \nN + 0(N°) , (1.3) 
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where r cusp (a s ) is the cusp anomalous dimension |21|. 

By definition, r cusp (o! s ) governs the scale dependence of Wilson lines with light-like 
cusps 1 22 ,23 1 and its relation to anomalous dimensions of large spin Wilson operators 
is by no means obvious. It can be understood |19| by invoking the physical picture of 
deep inelastic scattering at large N. In terms of the moments (1.1), large N corresponds 
to the region of x — > 1. For x — > 1 the final state in the deep inelastic scattering has 
a small invariant mass, s = Q 2 (l — x)/x <C Q 2 , and it consists of a collimated jet of 
energetic particles accompanying by soft gluon radiation. Interacting with soft gluons, 
the particles inside the jet acquire the eikonal phases given by Wilson line operators 
Pexp(i J °° dtp ■ A{pt)) evaluated along semi-infinite line in the direction of the particle 
momenta. In this way, for x — > 1, complicated QCD dynamics in deep inelastic scattering 
admits an effective description in terms of Wilson lines |24|. The relation (1.3) between 
anomalous dimensions and cusp singularities of light-like Wilson lines is just one of the 
application of this formalism. Another examples include the relation between light-like 
Wilson loops with on-shell scattering amplitudes, Sudakov form factors, gluon Regge 
trajectories etc (see Ref. |25| and references therein). 

At present, integrability of the dilatation operator in planar QCD has been verified to 
two loops in the SL(2; R) sector of maximal helicity operators |26|. In other sectors, the 
dilatation operator receives additional contribution that breaks integrability already to 
one loop. This contribution vanishes however for large values of the Lorentz spin N ^> 1 
thus suggesting that integrability in planar QCD gets restored to all loops in the leading 
large N limit |27|. Indeed, as was shown in Ref. |20|, the all-loop dilatation operator in 
QCD in the SL(2; R) sector can be mapped in the large N limit into a Hamiltonian of 
a classical Heisenberg SL(2;M.) spin chain. In this manner, the Wilson operators with 
large N are described by the so-called finite-gap solutions and the spectrum of anomalous 
dimension can be found through their semiclassical quantization. In particular, the 
relation (1.3) naturally appears as describing the ground state energy of the classical 
SL(2; M.) spin chain of an arbitrary length L and total spin N. 

The above mentioned integrability structures (those of the scattering amplitudes in 
the Regge limit and of the dilatation operator) are not specific to QCD. They are also 
present in generic four-dimensional gauge theories including supersymmetric Yang-Mills 
models with M = 1,2,4 supercharges. Supersymmetry enhances the phenomenon by 
extending integrability to a larger class of observables. In this context, the maximally 
supersymmetric M = 4 Yang-Mills theory is of a special interest with regards to the 
AdS/CFT correspondence |28|. The gauge/string duality hints that these structures 
should manifest themselves through hidden symmetries of the scattering amplitudes and 
of anomalous dimensions in dual gauge theories to all loops. 

2 Integrability of dilatation operator in QCD 

In this section, we review a hidden integrability of the dilatation operator in a generic 
four-dimensional Yang-Mills theory describing the coupling of gauge fields to fermions 
and scalars. Depending on the representation in which the latter fields are defined, we 
can distinguish two different types of the gauge theories: QCD and supersymmetric 
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extensions of Yang-Mills theory (SYM). 

In QCD, the gauge fields are coupled to quarks in the fundamental representation of 
the SU(N C ) gauge group. The quarks are described by four- component Dirac fermions 
ip and the gauge field strength F^ u = ^[D^,D U ] is determined in terms of the covariant 
derivatives = — igA a ^t a with generators t a in the fundamental representation of 
the SU(N C ) normalized conventionally as tr (t a t b ) = \b~ ab - In SYM theory, the gauge 
fields are coupled to fermions (gauginos) and scalars belonging to the adjoint represen- 
tation of the SU(N C ) group. The supersymmetric Yang-Mills theories with J\f = 1, 2 
and 4 supercharges are obtained from the Lagrangian of generic Yang-Mills theory by 
adjusting the number of gaugino and scalar species. The gauginos are described by the 
Weyl fermion X A which belongs to the fundamental representation of an internal SU{M) 
symmetry group with its complex conjugate = (\ A )*. The scalars are assembled 
into the antisymmetric tensor cj) AB = —<f) BA , with its complex conjugate ((p AB )* = 4>ab- 
As we explain below, integrability is not tied to supersymmetry and the phenomenon 
persists in the generic Yang-Mills theory for arbitrary A/", to two loop order at least. 



2.1 Light-ray operators 

Let us first consider renormalization of local gauge invariant operators in QCD. As 
the simplest example, we examine the following twist-two operator contributing to the 
moments of DIS structure function (1.1) 

(p\O N , L=2 {0)\p) = (p|^7+£>f"V(0)b> ■ (2.1) 

It is built from two quark fields and (N — 1) covariant derivatives D + = (n ■ D) pro- 
jected onto light-like vector n M = q^ — p^q 2 / (2pq) and 7 + = (n ■ 7) being the projected 
Dirac matrix. Discussing renormalization properties of Wilson operators like (2.1) it is 
convenient to switch from infinite set of local operators (2.1) parameterized by positive 
integer N to a single nonlocal light-ray operator 

(z — z ) N ^ 1 

<D(zi,z 2 ) = ^P(z 1 n) 1+ [nz 1 ,nz 2 Mz 2 n) = ]T ^ 7+ ^+"V] 1 ) N _ m +■■■ ( 2 - 2 ) 

N>1 ^ >' 

Here z\ and z 2 are scalar variables defining the position of quark fields on the light- 
cone and the gauge link [nzi,nz 2 ] = Pexp(ig f z 2 dtA + {nt)) is inserted to ensure gauge 
invariance of 0(21,22)- Also, ellipses in the right-hand side of (2.2) stand for terms 
involving total derivatives of the twist-two operators and, therefore, providing vanishing 
contribution to the forward matrix element (p|0(zi, z 2 )\p). 

We recall that local gauge invariant operators satisfy the evolution equation (1.2). 
The same is true for the light-ray operators (2.2) although the explicit form of the evo- 
lution equation is different due to nonlocal form of the light-ray operators. In particular, 
for the operators (2.2) the evolution equation takes the following form |29, 13,30,31 

^ + ^9 2 )^j @{zi,z 2 ) = -[m 2 (g 2 )-0}{ Zl ,z 2 ), (2.3) 
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with the evolution kernel H 2 to be specified below. The evolution equation (2.3) expresses 
the conformal Ward identity in QCD and the beta-function term takes into account 
conformal symmetry breaking contribution. The evolution operator H2 in the right-hand 
side of (2.3) defines a representation of the dilatation operator on the space spanned by 
nonlocal light-ray operators (2.2). In general, H2 has a matrix form as the light-ray 
operators with different partonic content could mix with each other. 

The evolution kernel EI 2 has a perturbative expansion in powers of the coupling con- 
stant and admits a representation in the form of an integral operator acting on light-cone 
coordinates Z\ and zi of 0(21,^2). To the lowest order in the coupling, its explicit form 
has been found in QCD in Ref. 1 30 1 and its generalization to Yang-Mills theories with 
an arbitrary number of supercharges has been derived in Ref. |32|. The corresponding 
expressions for H2 are given below in Eq. (2.11) . The main advantage of (2.3) compared 
with the conventional approach based on explicit diagonalization of the mixing matrix 
for local Wilson operators is that the problem of finding the spectrum of anomalous di- 
mensions can be mapped into spectral problem for one-dimensional quantum mechanical 
Hamiltonian H2. As we will see in a moment, the same happens in QCD for Wilson 
operators of high twist L > 3, in which case the corresponding evolution operator Ml in 
the sector of maximal helicity operators turns out to be equivalent for a Hamiltonian of 
Heisenberg SL(2; K) spin chain of length L. 

2.2 Light-cone formalism 

Discussing integrability of the dilatation operator in QCD and in SYM theories, it is con- 
venient to employ the "light-cone formalism" 1 33 - 35 1 . In this formalism one integrates 
out non-propagating components of fields and formulates the (super) Yang-Mills action 
in terms of "physical" degrees of freedom. Although the resulting action is not mani- 
festly covariant under the Poincare transformations, the main advantage of the light-cone 
formalism for SYM theories is that the M— extended supersymmetric algebra is closed 
off-shell for the propagating fields and there is no need to introduce auxiliary fields. This 
allows us to design a unifying light-cone superspace formulation of various M— extended 
SYM, including the M = 4 theory for which a covariant superspace formulation does not 
exist. 

In the light-cone formalism, one quantizes the Yang-Mills theory in a noncovariant, 
light-cone gauge (n ■ A) = A + (x) = 0. Introducing an auxiliary complimentary light-like 
vector n^, such that n 2 = and (n • n) = 1, we split three remaining components of the 
gauge field into longitudinal, A_(x), and two transverse holomorphic and antiholomor- 
phic components, A(x) and A(x), respectively, 

A_ = (n-A), A = j=(A l + iA 2 ), A = A* = j 3 (A 1 -iA 2 ). (2.4) 

In the similar manner, the fermion field ift(x) can be decomposed with a help of projectors 

^ = n+^ + n_^ = ^+ + ^-, (2.5) 

where the fermion field ip+ nas two nonzero components 

^ = ^(1-75)^+, qi = \{l + l^+- (2.6) 
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Then, one finds that the fields ip-(x) and A_(x) can be integrated out and the resulting 
action of the Yang-Mills theory is expressed in terms of "physical" fields: complex gauge 
field, A(x) and A(x), two components of fermion fields, q\(x) and q±(x), and, in the case of 
supersymmetric gauge theory, complex scalar fields 4>(x). When applied to the vacuum 
states, the fields (A, q±, <fi, q^, A) create massless particles of helicity (— 1, — 1,0, |,1), 
respectively. 

Taking the product of 'physical' fields and light-cone derivatives D + = d + , we can 
construct the set of local gauge invariant operators. Such operators define the represen- 
tation of the so-called collinear SL(2; K) subgroup of the conformal group and they are 
known in QCD literature as quasipartonic operators. A distinguished feature of these 
operators is that their twist equals the number of constituent physical fields |13|. In 
analogy with (2.2), we can replace an infinite number of Wilson operators of a given 
twist L with a few nonlocal light-ray operators 0(zi, . . . , Zl). The latter can be thought 
of as generating functions for the former. Due to different SU(N C ) representation of 
fermions (fundamental in QCD and adjoint in SYM), the definition of such operators is 
slightly different in the two theories. 

In QCD, in the simplest case of twist two, we can distinguish four different light-ray 
operators (plus complex conjugated operators) 

0^( Zl , z 2 ) = q^nz x )q^nz 2 ) , 0$( Zl , z 2 ) = tr [d + A(n Zl )d + A(nz 2 )] , 

<□>«(*!, z 2 ) = q^nz^nz*) , O gg \ Zl , z 2 ) = tr [d+A(nzi)d + A(nz 2 )} , (2.7) 

where the subscript (qq and gg) indicates particle content of the operator and the su- 
perscript defines the total helicity. In this basis, the operator (2.2) is given by a linear 
combination of O q °q(zi,z 2 ) and complex conjugated operator. The operators Oqq and 
Ogg have the same quantum numbers and mix under renormalization. At the same time, 
the operators 0$ and OgJ carry different helicity and have an autonomous scale depen- 
dence. In what follows we shall refer to them as maximal helicity operators. The reason 
why we distinguish such operators is that the one-loop dilatation operator in QCD is 
integrable in the sector of maximal helicity operators only. 

For higher twist L > 3 we can define three different types of maximal helicity opera- 
tors in QCD: 

®2 2) (2i, *a, zs) = e ijk q\{ Zl n)q{{z 2 n)q^{z 3 n) , (2.8) 
Ofc$(*i> • • • > z l) = qi(n Zl )d + A(nz 2 ) . . . d + A(nz L ^)q^nz L ) , (2.9) 
0( g L } g (z 1 ,...,z L )=tr[d + A(nz 1 )...d + A(nz L )) , (2.10) 

to which we shall refer as baryonic (L = 3) operators, mixed quark-gluon operators and 
gluon operators, respectively. We remind that since quark fields belong to the funda- 
mental representation of the SU(N C ) group, the length of the operator (2.8) ought to be 
N c = 3. At the same time, gluon fields are in the adjoint representation and the single 
trace operator (2.10) is well-defined for arbitrary N c and twist L. The same applies 
to the mixed quark-gluon operators (2.9). The operators (2.8) and (2.9) have a direct 
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phenomenological significance: their matrix elements determine the distribution ampli- 
tude of the delta-isobar |36| and higher twist contribution to spin structure functions, 
respectively. 



2.3 Evolution kernels 



The light-ray operators (2.7) - (2.10) satisfy the evolution equation (2.3). Let us first 
examine twist-two quark operators and defined in (2.7). The operator 
can mix with the gluon operator O g °g. To simplify the situation, we can suppress the 
mixing by choosing the two quark fields inside to have different flavor. To one- loop 
order, the evolution kernel receives the contribution from one-gluon exchange between 
two quark fields and from self-energy corrections. The latter one is the same for the two 
operators while the former one is different 

,2 



M ) = ^%[tf 12 +V 12 +2 T J . (2.11) 
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Here C F = t a t a = (iV c 2 - 1)/(2JV„) is the quadratic Casimir of the SU(N C ) in the fun- 
damental representation, ^ q = 1 is one-loop anomalous dimension of quark field in the 
axial gauge A + = and Hi 2 and Vi 2 are integral operators 



da _ 
— a 
a 



Zi, z 2 ) - U(az 1 + az 2 , z 2 ) - ( u>{z 1 , az x + az 2 ) 

1 pai 



{H 12 -0]( Zl ,z 2 ) = 

[Vi 2 ■ 0)(zi, z 2 ) = / da 1 / da 2 <0(aiz 1 + aiz 2 ,a 2 z 2 + a^) , (2.12) 
Jo Jo 



where cii = 1 — ctj. These operators have a transparent physical interpretation: they 
displace two particles along the light-cone in the direction of each other. 

To find the spectrum of anomalous dimensions of twist-two quark operators generated 
by light-ray operators (2.7), we have to diagonalize the operators W$ and M. q °J . This 
can be done with a help of conformal symmetry. We recall that the conformal symmetry 
is broken in QCD at loop level. However the dilatation operator receives conformal 
symmetry breaking contribution only starting from two loops and, as a consequence, 
the one- loop evolution kernels in QCD have to respect conformal symmetry of QCD 
Lagrangian. For nonlocal light-ray operators built from fields X(nz), the full 50(2,4) 
conformal symmetry reduces to its collinear SL(2; R) subgroup acting on one-dimensional 
light-cone coordinates of fields |14,15| 



az + b 

z ->■ 



X(zn) -> (cz + d)~ 2j X ( ^^-n) (2.13) 

\cz + d J 



CZ + G? ' 

with ad — be = 1. The generators of these transformations are 

L- = -d z , L+ = 2jz + z 2 d z , L°=j + zd z . (2.14) 
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Here j is the conformal weight of the field. For 'physical' components of fermions, ip + , 
it equals j q = 1, for transverse components of gauge field, d + A and d + A, it is j g = 3/2 
and for scalars j s = 1/2. 

In application to light-ray quark operators, O qq (zx, z 2 ) and O qq \zi, z 2 ) } the conformal 
symmetry dictates that the one-loop evolution kernels (2.11) have to commute with the 
two particle conformal spin + L 2 (with a = — , +, 0). As a consequence, H^ -0 ' 1 ^ is a 
function of the corresponding two-particle Casimir operator 

L 2 12 = (L a 1+ L a 2 ) 2 = J 12 {J 12 -l). (2.15) 

Q = + ,-,0 

To find the explicit form of this dependence, it suffices to examine the action of the 
two operators, and Lf 2 , on the same test function [z\ — z 2 ) n , which is just the 
lowest weight in the tensor product of two SL(2; IR) representations carrying the spin 
J 12 = n + 2. Replacing 0(zi, z 2 ) ->■ {z x - z 2 ) Jl2 ~ 2 in (2.12) we find 

H 12 = 2 [^(Ji 2 ) - ^(2)] , V 12 = 1/(J 12 (J 12 - 1)) , (2.16) 

where ip(x) = dlnT(x)/dx is Euler psi-function. Together with (2.11) these relations 
determine the spectrum of anomalous dimensions of twist-two quark operators. 



2.4 Relation to Heisenberg SX(2;IR) spin chain 

As the first sign of integrability, we notice that H\ 2 coincides with the known expression 
for two-particle Hamiltonian of Heisenberg spin chain |37,38| 

H L = H 12 + ... + H L1 , F M+1 = VOW) - V(2j) , (2.17) 

where the spin operators are identified as SX(2;R) conformal generators (2.14). As fol- 
lows from (2.11), the one-loop dilatation operator depends on H i2 and, therefore, 
it is mapped into Heisenberg 5L(2;R) spin chain of length 2. At the same time, the 
dilatation operator M qq receives the additional contribution Vi 2 . It preserves the confor- 
mal symmetry but breaks integrability. Notice that V\ 2 vanishes for large values of the 
conformal spin J\ 2 ^> 1 so that the two evolution kernels, Mqq and M. qg \ have the same 
asymptotic behavior at large J\ 2 . This suggests that for the operator W qq integrability 
is restored in the limit of large conformal spin only. 

For twist-two operators, the anomalous dimensions are uniquely determined by their 
conformal spin. To appreciate the power of integrabilty, we have to consider Wilson 
operators of high twist L > 3. For example, for the maximal helicity baryonic operators 
(2.8) the one-loop dilatation operator has the form |10| 



H (3/2) = ^ 

qqq 2vr 



1/N C )(H 12 + H 23 + H 31 ) + ^C F 



(2.18) 



with 7V C = 3 and Hi 2 given by (2.16). Comparing this relation with (2.17) we recognize 

(3/2) 

that Wqq^ can be mapped into Heisenberg spin chain of length L — 3. The spin at each 
site j = 1 is determined by the conformal spin of quark field. 
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For gluon operators of the maximal helicity (2.10) the dilatation operator receives 
contribution from self-energy corrections to gluon fields and from one-gluon exchange 
between any pair of gluons. The latter produces both planar and nonplanar corrections 
(for L > 3). In the planar limit, the one-loop dilatation operator has the following 
form 1 39 

^ 9 L) , = 9 -^(H l2 + ... + H L1 ) } (2.19) 

where two-particle kernel i^+i acts locally on light-cone coordinates of gluons with 
indices i and % + 1. The conformal symmetry implies that i/^i+i is a function of the 
conformal spin of two gluons Jij+i- Quite remarkably, the dependence of H iti+ i on 
J^i+i has the same form as in (2.17). As a consequence, the one-loop planar dilatation 
operator for maximal helicity gluon operator (2.10) coincides with the Hamiltonian of 
the Heisenberg SL(2; R) spin chain. The length of the spin chain equals the twist of the 
operator L and the spin in each site j = 3/2 coincides with the conformal spin of the 
gluon field. 

For mixed quark-gluon operators of the maximal helicity (2.9), the quark fields can 
interact in the planar limit with the adjacent gluon fields only while quark-quark inter- 
action is suppressed in this limit. As a consequence, the one- loop dilatation operator has 
the following form in the planar limit 

M i L g :.ll = ^(Ui2 + H 23 + ... + H L _ 1>L + U L _ 1>L ) . (2.20) 

Here H i i+ i describes the interaction of two gluons with aligned helicities and it is the 
same as in (2.17). The kernels Ui 2 and Ul-i,l describes quark-gluon interaction and 
their explicit form can be found in Ref. |40,41|. Notice that the operator M.q^..^\ has the 
form of a Hamiltonian of open spin chain of length L. The spin in sites 1 and L coincides 
with the conformal spin of quark j q = 1 and the spin in all remaining sites is given by 
gluon conformal spin j g = 3/2. As was shown in Ref. |40,41|, the open spin chain (2.20) 
is integrable. 



2.5 Exact solution 



Integrability of the one-loop dilatation operator allows us to find the exact spectrum of 
anomalous dimensions with a help of the Bethe Ansatz |8-10 

9 2 N C 4 

lN,L = -^-j-ZN,L + 0{g ) , 



£ 



N,L 



2j 



N 

y 

f^ u l+f du Q(u-tj 



■ d Q(u + ij) 
i— In 



(2.21) 



u=0 



Here j is the conformal spin in each site (j = 1 for quark operators and j = 3/2 for 
gluon operators), Uk are Bethe roots and Q{u) is a polynomial of degree N of the form 



Q(u) 



N 



Uj ). 



(2.22) 
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The function Q(u) defined in this way has the meaning of the eigenvalue of the Baxter 
operator for the SL(2, R) magnet |4,42|. It satisfies the finite-difference Baxter equation 

t L {u)Q{u) = {u + i]) L Q{u +i) + (u- i]) L Q{u - i) , (2.23) 

where ti{u) is the transfer matrix of the spin chain 

t L (u) = 2u L + q 2 u L - 2 + ... + q L (2.24) 

and q 2 , . . . ,qL are the conserved charges. 

The Baxter equation (2.23) alone does not specify Q(u) uniquely and it has to be 
supplemented by additional condition for analytical properties of Q(u). For the SL(2; R) 
spin chains describing the anomalous dimensions, Q(u) has to be a polynomial in the 
spectral parameter. Being combined with the Baxter equation (2.23), this condition 
determines Q(u) up to an overall normalization and, as a consequence, allows us to 
establish the quantization conditions for the q— charges and to compute the exact energy 

£-N,L- 

Solving the Baxter equation (2.23) for N = 0,1, . . . one finds the eigenspectrum of the 
Hamiltonian and, as a consequence, determines the exact spectrum of the anomalous 
dimensions of the maximal helicity baryon operators (for j = 1 and L = 3) and of 
maximal helicity gluon operators (for j = 3/2 and L > 2). The spectrum obtained in 
this way exhibits remarkable regularity: almost all eigenvalues are double degenerate and 
for large N they belong to the set of trajectories |4,43|. Both properties are ultimately 
related to integrability of the dilatation operators and can be served to test integrability 
at high loops. 

For the SL(2; R) spin chains under consideration, the Baxter equation approach and 
conventional Bethe Ansatz are equivalent. Indeed, substituting (2.22) into the Baxter 
equation (2.23), one finds that the roots Uj satisfy the conventional SU{2) Bethe equa- 
tions for spin (—j). The fact that the spin is negative leads to a number of important 
differences as compared to "compact" SU{2) magnets. In particular, the Bethe roots 
take real values only and the number of solutions is infinite |4,43|. 



2.6 Semiclassical limit 

The Baxter operator approach becomes advantageous when one studies the properties 
of anomalous dimensions at large spin N and/or twist L. The reason for this is that 
the Baxter equation (2.23) takes the form of discretized Schrodinger equation. After 
rescaling of the spectral parameter, u — > (N + Lj)x, we can seek for solution to (2.23) 
in the WKB form 44,45,43 

Q(Nx) = exp (is(x)j , h=l/(N + Lj), (2.25) 

where the action function S(x) admits an expansion in powers of h. Substitution of 
(2.25) into the Baxter equation (2.23) yields the equation for S(x) which can be solved 
as a series in h. To leading order we have 

S(x)= dxp(x) + 0(h), (2.26) 

Jx 



11 



where the momentum p(x) is defined on the spectral curve ("equal energy" condition) 
of the classical SX(2;R) magnet y(x) = 2x L sinhp(x) with |46 

T L : y 2 = {t L {x)f - Ax 2L . (2.27) 

The classical dynamics on this spectral curve has been studied in detail in Refs. |47,27|. 
Using (2.25) we can compute the asymptotic behavior of the energy as |43,48 

L 

4 a J = 2 In 2 + ]T + iS n ) + - iS n ) - ${2j)\ + ..., (2.28) 

n=l 

where ellipses denote terms subleading at large (N+jL). Here 5 n are roots of the transfer 
matrix defined in (2.24), ti(S n ) = 0. They depend on the conserved charges q2,---,QL 
whose values satisfy the WKB quantization conditions 

* dxp{x) = 27i7i(4 + §), (for k = 1,...,L-1). (2.29) 

Here integration goes over the cycles on the complex curve (2.27) encircling intervals 
on the real axis satisfying y 2 (x) > and integers £k enumerate the quantized values of 
the charges q^, . . . , qi and the energy £m,l = £n,l(£i, ■ ■ ■ ,^1-2)- For large spin iV and 
twist L, the minimal energy £jtp L = min£ fc £n,l has the following scaling behavior |48 

4% = f(p) In N + 0(N°), p =J- = £ X ed, N, L » 1 (2.30) 

where f(p) is the so-called generalized scaling function. Detailed analysis of the relations 
(2.28) and (2.29) can be found in Refs. 143,10,39,27,481- For recent development in the 
generalized scaling function in M = 4 SYM see review Ref. |49|. 

So far we have discussed the exact solution for the one-loop anomalous dimensions 
of quark and gluon maximal helicity operators. For the anomalous dimensions of mixed 
quark-gluon operators (2.9), similar analysis of the spin chain (2.20) can be carried out 
using Bethe Ansatz for open SL(2;M) spin chains |39,40|. 

2.7 Integrability of dilatation operators in SYM theories 

In this subsection, we extend consideration to supersymmetric Yang-Mills theories. Dis- 
cussing integrability of dilatation operator in these theories, it is convenient to employ 
supersymmetric version of light-cone formalism due to Mandelstam |35| and Brink et 
al. |34|. In this formalism, all symmetries of SYM theory become manifest and cal- 
culations can be performed in a unified manner for different numbers of supercharges 
M = 0,1,2,4. The maximally-supersymmetric M = 4 SYM theory is a finite, four- 
dimensional conformal field theory 1 35, 34, 50, 51 1, while the M = theory corresponds 
to pure gluodynamics. 

Defining a SYM theory on the light-cone, one starts with the component form of the 
action, fixes the light-cone gauge A + (x) = 0, decompose all propagating, "physical" fields 
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into definite helicity components. In the case of M = 4 SYM, they include helicity (±1) 
fields, A(x) and A(x), built from two-dimensional transverse components of the gauge 
field, complex scalar fields <p AB of helicity and helicity ±1/2 components of Majorana- 
Weyl fermions, \ A and Xa, all in the adjoint representation of the SU(N C ) gauge group. 
An important property of the light-cone formalism, which makes it advantageous over 
the covariant one, is that the latter fields have only one non-vanishing component. As 
a consequence, one can describe helicity (±1/2) fermions by Grassmann- valued complex 
fields without any Lorentz index. Introducing four fermionic coordinates 9 A (with A = 
1,...,4) possessing the helicity (— |) and their conjugates 9 a with helicity \, we can 
assemble the above fields into a single, complex chiral M = 4 superfield |34 

${x,9 A ) = d+ 1 A{x) + 9 A d+ 1 \ A {x) + ^9 A 9 B AB {x) 

- ^e A BCD9 A 9 B 9 c \ D {x)-^e A BCD9 A 9 B 9 c 9 D d + A{x). (2.31) 

It embraces all particle helicities, from —1 to 1 with half-integer step, and, therefore, 
Q(x,9 A ) describes a CPT self-conjugate supermultiplet. 

Gauge theories on the light-cone with less or no supersymmetry can be deduced from 
the maximally supersymmetric M = 4 theory by removing "unwanted" physical fields. 
In the superfield formulation this amounts to a truncation of the Af = 4 superfield, or 
equivalently, reduction of the number of fermionic directions in the superspace |52|. For 
instance, to get the Af = 1 superfields one removes three odd coordinates 9 2 = 9 3 = 
9 A = 0, whereas for M = all 9's in (2.31) have to be set to zero. Notice that under this 
procedure the truncated M = 2, M = 1 and M = theories involve only half of the fields 
described by the Af— extended SYM theory and the other half of the needed particle 
content arises from the complex conjugated superfields $ = $*. Explicit expressions 
for the action of the SYM theory in terms of the light-cone superfields can be found in 
Ref. [32]. 

In a close analogy with (2.10), we can introduce multiparticle single-trace operators 
built from light-cone superfields 

0(Z h ... } Z L )=tr {$(Z 1 )$(Z 2 ) • • • $(Z L )} , (2.32) 

where $(Z) = § a (Z)t a is a matrix (SU(N C )) valued superfield and Z = (x,9 A ) denotes 
its position in the superspace with four even coordinates, x M , and Af odd coordinates, 
9 A with A = 1, . . . ,Af. In addition, we choose all superfields to be located along the 
light-cone direction in the four- dimensional Minkowski space defined by the light-like 
vector n M (with n 2 = 0), so that n ■ A = A + = 0. Similarly to the QCD case, the posi- 
tions of the superfields on the light-cone are parameterized by real numbers x M = zn^, 
$(Zfc) = $(zfcn, 9 A ) . The single-trace operators (2.32) represent a natural generalization 
of nonlocal light-ray operators in QCD, cf. Eq. (2.10). To obtain the latter it is sufficient 
to expand 0(Z 1; . . . , Z{) in powers of odd variables 9 Al . . . 9 Al . As in QCD, nonlocal 
operators (2.32) serve as generating functions for Wilson operators with the maximal 
Lorentz spin and minimal twist equal to the number of constituent fields L. Such oper- 
ators define a representation of the SL(2\Af) subgroup of the full superconformal group. 
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Examining light-ray operators (2.32) in SYM theories with different number of su- 
percharges, we find that M = 4 case is special. In M = 4 SYM theory there is only one 
independent chiral superfield &(Z) and, as a consequence, the operators (2.32) generate 
all Wilson operators of twist— L built from L fundamental fields. For J\f < 2, the su- 
perfield $(Z) and its conjugate &(Z) are independent of each other and, in addition to 
the operators in (2.32), one can introduce "mixed" operators built from both superfields. 
This means that in the J\f — 0, 1 and 2 SYM theories, the operators (2.32) only generate 
a certain subset of the existing Wilson operators in the SL(2\Af) subsector. 

The light-ray operators (2.32) play a special role as far as integrability is concerned. 
Namely, as was shown in Refs. |32|, the one- loop dilatation operator acting on the space 
of single-trace operators (2.32) can be mapped in the multicolor limit into a Hamiltonian 
of a completely integrable Heisenberg SL(2\M) spin chain. As before, the length of the 
spin chain coincides with the number of superfields in (2.32) and spin operators are 
generators of a collinear SL{2\M) subgroup of the full superconformal group |32|. 

We recall that in SYM theories with J\f < 2 supercharges the operators (2.32) 
only generate a subsector of Wilson operators of twist L. To describe the remaining 
operators, one has to consider single-trace operators built from both superfields, like 
tr {$(Zi)$(Z2) • • • $(Z^)}. For such operators, the one-loop dilatation operator involves 
the additional term describing the exchange interaction between superfields on the light- 
cone — > $$. It breaks integrability symmetry and generates a mass gap in the 
spectrum of the anomalous dimensions |10|. At the same time, for large values of the 
superconformal spin the exchange interaction vanishes and integrability gets restored in 
the leading large spin asymptotics of the anomalous dimensions. 

2.8 Integrability in QCD and SYM beyond one loop 

It is well-known that the conformal symmetry is broken in QCD and SYM theories with 
M < 4 supercharges while in the maximally supersymmetric M = 4 model it survives 
on the quantum level. However the conformal anomaly modifies anomalous dimensions 
starting from two loops only and, therefore, the one-loop dilatation operator inherits the 
conformal symmetry of the classical theory |11,53|. 

Starting from two-loop order, the dilatation operator in the SL(2) sector acquires 
several new features. First, it receives conformal symmetry breaking corrections arising 
both due to a nonzero beta-function and a subtle symmetry-violating effect induced by 
the regularization procedure |54|. Second, the form of the dilatation operator starts to 
depend on the representation of the fermion fields, i.e., fundamental SU(3) in QCD and 
adjoint SU(N C ) in SYM theories. The difference between the two is that it is only in the 
latter case that one can select planar diagrams by going over to the multi-color limit, 
while in the former case the large- N c counting is inapplicable and the two-loop dilata- 
tion operator receives equally important contributions from both planar and nonplanar 
Feynman graphs. Thus, by studying the two-loop dilatation operator in the SL(2) sec- 
tor we can identify what intrinsic properties of gauge theories (conformal symmetry, 
supersymmetry and/or planar limit) are responsible for the existence of the integrability 
phenomenon per se. 

For an all-loop dilatation operator HI(A), depending on 't Hooft coupling constant 



14 



A = g 2 N c /(8ir 2 ) and acting on a Wilson operator built from L constituent fields and 
an arbitrary number of covariant derivatives, integrability would require, in general, the 
existence of L conserved charges. Two of the charges — the light-cone component of the 
total momentum of L fields and the scaling dimension of the operator — follow immedi- 
ately from Lorentz covariance of the gauge theory. However, the identification of the 
remaining charges ^(A) with k = 3, . . . , L is an extremely nontrivial task. The eigen- 
values of the charges define the complete set of quantum numbers parameterizing the 
eigenspectrum of the dilatation operator. Integrability imposes a nontrivial analytical 
structure of anomalous dimensions of Wilson operators and implies the double degen- 
eracy of eigenvalues with the opposite parity |4,55,10|. At the same time, breaking 
of integrability leads to lifting of the degeneracy in the eigenspectrum of the one-loop 
dilatation operator. 

Explicit two-loop calculation of the anomalous dimensions of the aforementioned 
aligned-helicity fermionic operators in all SYM theories showed that the same relation 
between integrability and degeneracy of the eigenstates holds true to two loops. Namely, 
as was found in Refs. |26|, the desired pairing of eigenvalues occurs for three-gaugino 
operators in SYM theories with Af = 1, 2 supercharges and the SU(N C ) gauge group. 

The two-loops dilatation operator in SYM theories receives conformal symmetry 
breaking contribution and, in addition, it depends on the number of supercharges Af '. 
The latter dependence comes about through the contribution of 2 (Af— 1) real scalars and 
Af gaugino fields propagating inside loops. Both contributions to two-loop dilatation op- 
erator can be factored out (modulo an additive normalization factor) into a multiplicative 
c-number. This property makes the eigenspectrum of the two-loop dilatation operator 
alike in all gauge theories including the Af = 4 SYM in which case the dilatation operator 
is believed to be integrable to all loops |56|. Summarizing the results of two- loop cal- 
culations of the anomalous dimension in QCD and in SYM theories, integrability of the 
dilatation operator only requires the planar limit but it is sensitive neither to conformal 
symmetry, nor to supersymmetry |26|. For recent discussion of integrability in relation to 
non-planar corrections to the anomalous dimensions in Af = 4 SYM see review Ref. |57|. 

In this section, discussing the properties of anomalous dimensions we restricted our- 
selves to the SL(2) sector. There have been several developments that we cannot address 
here in detail. In particular, an important observation was made in Ref. |58|, where it 
was shown that the diagonal part of one-loop QCD evolution kernels governing the scale 
dependence of Wilson operators of arbitrary twist, can be written in a Hamiltonian form 
in terms of quadratic Casimir operators of the full conformal 5*0(2, 4) group. This obser- 
vation was used in Ref. |59| to work out the non-diagonal parts of the evolution kernels 
for generic twist-four operators. 

3 Integrability in high energy scattering 

In the previous section, we described how integrability emerges in the problem of finding 
the dependence of the structure functions F(x, Q 2 ) on the hard scale Q 2 . In this section, 
we explain that yet another integrability symmetry arises in the high-energy limit. 
In application to the structure function F(x, Q 2 ) this limit corresponds to x — > for 
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fixed Q 2 . At small x, the invariant energy s = Q 2 (l — x)/x of colliding virtual photon and 
hadron becomes large and the structure function is expected to have Regge-like scaling 
behavior F(x, Q 2 ) ~ In terms of moments (1.1), this corresponds to appearance 

of the Regge pole at N = oj 

F N {q 2 ) = f dxx N - l F{x,q 2 ) ~ - J— . (3.1) 

JO 1\ — LO 

It is well-known |5| that perturbative corrections to F(x, q 2 ) are enhanced at small x by 
large logarithms ~ (a s hi(l/x)) p . This raised the hope that the Regge behavior (3.1) can 
be derived in QCD from resummation of such corrections to all loops. Going to moments, 
the expansion over (a s lns) p is traded for the expansion of F^(q 2 ) over (a s /N) p . 

3.1 Evolution equation 

Careful study of asymptotic behavior of Feynman diagrams describing interaction be- 
tween virtual photon and hadron shows that the dominant contribution to F(x, q 2 ) only 
comes t— channel exchange of particles of maximal spin, i.e. gluons (see Fig. 2). More- 
over, in the center-of-mass frame of 7* (5) and h(p), due to hierarchy of the scales, s > Q 2 , 
interaction takes place in the two-dimensional plane orthogonal to the plane defined by 
the momenta of scattered particles, and q^. This implies that in generic Yang-Mills 
theory the leading high-energy asymptotic behavior of the scattering amplitudes is driven 
by t— channel exchange of an arbitrary number of gluons. In the so-called generalized 
leading logarithmic approximation, their contribution to the moments (3.1) takes the 
form 

F N (q 2 ) = J2 fW f[d 2 k'] <f> r ({k})T L ({k},{k'};N)$ h ({k'}) , (3.2) 

L>2 J J 

where integration goes over two-dimensional momenta of L gluons propagating in the 
t— channel, [d 2 k] = Y\i d 2 ki and similarly for [gP/c']. Here, the wave functions $ 7 *({/c}) 
and &h({k'}) describe the coupling of L gluons to virtual photon and hadron, respectively. 
Also, Ti({fc}, {&'}; N) describes elastic scattering of L gluons in the t— channel (see 
Fig. 2) and is the main object of our consideration. 

It is convenient to rewrite (3.2) as the following matrix element 

F N (q 2 ) = J2(®r\TL(N)\®h), (3.3) 

L>2 

where the minimum number of two gluons, L = 2, is required in order to get a colorless 
exchange. The transition operator Tl(N) describes the elastic scattering of L gluons. In 
the generalized leading logarithmic approximation, the Feynman diagrams contribution 
to Tl(N) have ladder structure as shown in Fig 2. They can be resummed leading to 
the following Bethe-Salpeter equation |6,7 

NT L (N) = Tf + ^U L T L (N) , (3.4) 
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Figure 2: The Feynman diagrams contributing to the deep inelastic scatter- 
ing in the generalized leading logarithmic approximation. Wavy lines denote 
(reggeized) gluons. They couple to virtual photons through a quark loop. 

where corresponds to the free propagation of L gluons in the t— channel and the 
evolution operator Ml describes their pair- wise interaction. The operator M L acts both 
on two-dimensional momenta and on colors of L gluons and has the following two-particle 
form 

U L = J2 HiiWl ( 3 - 5 ) 

l<i<j<L 

Each term in this sum is given by the product of the color factor involving color charges 
of two gluons and two-particle kernel if^ acting locally on the tranverse momenta of 
gluons with indices i and j. The kernel is known as BFKL operator |5| and it is 
denned below in (3.19). 

Combining together (3.4) and (3.3) we obtain the following expression for F]y(q 2 ) 

Ml 2 ) = J> 7 *|(iV- ^H L ) -1 7f \9 h ) . (3.6) 

L>2 

We observe that F N (q 2 ) has (Regge) singularities in N which are determined by the 
eigenspectrum of the operator M n , the so-called BKP equation |6,7|, 

H L ^L,{ g #i, ...,k L ) = EL^L^yih, ...,k L ). (3.7) 

The solutions to (3.7) define color singlet compound states of L gluons and we introduced 
{q} to denote the set of quantum numbers parameterizing all solutions. Having solved 
Schrodinger like equation (3.7), we can compute the moments of the structure function 

as |4 

^v(<7 2 ) = |^l,m) ~Vl, W ■ (3.8) 

L>2 { q } 

Here the impact factor = (^-r*\^L ; {q})(^L,{q}\T^\^h} measures the projection 

of the eigenstates onto the wave functions of scattered particles. The double sum in 
(3.8) runs over the possible number of gluons L > 2 and over all eigenstates of the 
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BKP Hamiltonian (3.7) parameterized by the conserved charges q. We observe that this 
relation has an expected Regge form (3.1). Moreover, the leading Regge behavior of the 
structure function is controlled by right-most singularity of F^{q 2 ) in complex N plane. 
According to (3.8), it corresponds to the maximal value of the 'energy' El,{ q }- 



3.2 Conformal SL(2; C) symmetry 

We recall that ki in the BKP equation (3.7) describe two-dimensional transverse momenta 
of ith gluon and the relation (3.7) can be interpreted as two-dimensional Schrodinger 
equation for n particles carrying SU(N C ) color charges. 

As was found in |60,2,3|, the BKP equation (3.7) becomes integrable in the multi- 
color limit. In this limit, the relevant ladder Feynman diagrams contributing to F^(q 2 ) 
have the topology of a cylinder and, as a consequence, the evolution operator Ml reduces 
to the sum of terms corresponding to pairwise nearest-neighbor BFKL interactions: 

1 L 

M L = -Y,H k , k+l + 0{l/N 2 c ) (3.9) 

k=l 

with periodic boundary conditions Hl,l+i — Hl,i- Notice that this relation is exact for 
L = 2. 

The BFKL operator H ktk+ i has a number of remarkable properties which allow us to 
solve the Schrodinger equation (3.7) exactly |61,60|. To elucidate these properties it is 
convenient to switch from two-dimensional momenta ki to two-dimensional coordinates bi 
via Fourier transform and, then, introduce complex holomorphic and the antiholomorphic 
coordinates 

— * —> 

ki i-> bi = {xi,yi} !->■ (zi = Xi + iyi, z { = x { - iyi) . (3.10) 

Quite remarkably, H\ 2 is invariant under the conformal SL(2; C) transformations of the 
gluon coordinates on the plane 1 6 1,60 

CLZ k ~ b . . 

z k -)■ — , (ad-bc=l), (3.11) 

cz k + a 

and similarly for antiholomorphic coordinates z k . The generators of these transforma- 
tions are 

L k - = —d Zk , L k fi = z k d Zk , L fc + = z\d Zh , (3.12) 

and the corresponding antiholomorphic generators L k _, L k} o and L kj+ are given by sim- 
ilar expressions with z k replaced by z k , with k = 1,2 enumerating particles. Then, i7 12 
commutes with all two-particle generators 

[#12, L 1>a + L 2>a ] = [H 12 , L lA + Z 2)0 ] = (3.13) 

with a = +, — , 0. This implies that, firstly, if 12 only depends on the two-particle Casimir 
operators of the SL(2, C) group 

L i2 = - z 2 fd Zl d Z2 , L\ 2 = -(zi - z 2 fd Zl d l2 , (3.14) 
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and, secondly, the eigenstates of Hyi have to diagonalize the Casimir operators 



I&a ntV = h{h - l)* nji/ , L%& n>v = h{h-\)V n j / . (3.15) 
Here a pair of complex conformal spins is introduced 

1 -\~ TI — 1 Ti 

h = — \-iu, h = — h iv (3.16) 

with a non-negative integer n and real v that specify the irreducible (principal series) 
representation of the SL(2, C) group to which ^ n%v belongs. The solutions to Eqs. (3.15) 
are |61 

\ (l+n)/2+w , - \ (l-n)/2+«/ 
%12 \ I %12 



*n,u{K 62) = — ^~ ~ ^" , (3-17) 

\ z 10^20/ \ Z 10^20/ 

where Zjk = Zj — Zk and bo = (za,Zo) is the collective coordinate, reflecting the invariance 
of H12 under translations. The corresponding eigenvalue of H\ 2 reads |5,61 

E n „ = 2^(1) - i> +iv)-^ - iv) . (3.18) 



Its maximal value, ma.xE UjU = 4 In 2, corresponds to n = v = 0, or equivalently h — h = 
1/2. It defines the position of the right-most singularity u = 4 In 2 a s N c / n in (3.1 ) known 
as the BFKL pomeron |5|. The relations (3.17) and (3.18) define the exact solution to 
the Schrodinger equation (3.7) for n = 2, that is for the color-singlet compound state 
built from two reggeized gluons. 



3.3 Heisenberg SL(2; C) spin chain 

Using (3.18) one can reconstruct the operator form of the BFKL kernel H12 on the 
representation space of the principal series of the SL(2, C) group 

H 12 = l - [H(J 12 ) + H{J 12 )] , H(J) = 2^(1) - ip(J) - V(l - J) , (3.19) 

where, as before, the two-particle spins are defined as L\ 2 = J\ 2 {Ji2 — 1) and L\ 2 = 
Jyi{Jyi — !)■ Notice that we already encountered the similar Hamiltonian in Sect. 2 (see 
Eq. (2.17)) and found that it gives rise to integrability for the dilatation operator. 

Most remarkably, the Hamiltonian (3.9) has the same hidden integrability as the 
dilatation operator (2.17) and it coincides in fact with the Hamiltonian of the SL(2,C) 
Heisenberg magnet |2,3|. The important difference between the two operators is that 
they are defined on the different space of functions: the operator (2.17) acts on the 
nonlocal light-ray operators (2.2) which are polynomials in the light-cone coordinates 
while the eigenfunctions of the operator (3.19) are single- valued functions on the two- 
dimensional plane, Eq. (3.17). This leads to a dramatic change in the properties of the 
two evolution kernels. 

The number of sites in the Heisenberg SL(2, C) spin chain (3.9) equals the number of 
particles and the corresponding spin operators are identified as six generators, L^, L° and 
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L^,L°, of the SX(2,C) group. It possesses a large-enough set of mutually commuting 
conserved charges q n and q n (n = 2, L) such that q n = q\ and [H L , q n ] = \M. L , q n ] = 0. 
The charges q n are polynomials of degree n in the holomorphic spin operators. They 
have the following form 1 2 , 3 

q n = Z jlj2 Z j2j3--- Z jnjlPjlPj2---Pjn (3.20) 

l<h<32<-<j n <L 

with zjk = Zj — Zk and pj = id Zj . The "lowest" charge q2 is related to the total spin of 
the system h. For the principal series of the SL(2, C) it takes the following values 

q 2 = -h(h-l), h = ^^ + iv h , (3.21) 

with rih integer and Vh real. The eigenvalues of the integrals of motion, 52, form 
the complete set of quantum numbers parameterizing the L— gluon states (3.7). 

Identification of (3.9) as the Hamiltonian of the SL(2,C) Heisenberg magnet allows 
us to map the L— gluon states into the eigenstates of this lattice model. In spite of the 
fact that the Heisenberg SX(2,C) magnet represents a generalization of the SL(2,M.) 
spin chain, finding its exact solution is a much more complicated task. The principal 
difficulty is that, in distinction with SL(2, K) magnet, the quantum space of the SL(2, C) 
magnet does not possess a highest weight - the so-called "pseudo-vacuum state" - and, 
as a consequence, conventional methods like the Algebraic Bethe Ansatz method |16| are 
not applicable. The eigenproblem (3.9) has been solved exactly in Refs. |62-64| using 
the method of the Baxter Q— operator 1 65, 45 ,66, 44 1 which does not rely on the existence 
of a highest weight. In this approach, it becomes possible to establish the quantization 
conditions for the integrals of motion (fe, ...,c?z, and to obtain an explicit form for the 
dependence of the energy El on the integrals of motion. 

In this manner, the spectrum of the L— gluon state has been calculated for L > 3 
particles: For L = 3 few low-lying states have been found in |67,68| and the complete 
spectrum of states for 3 < L < 8 was determined in [63 , 64 1 (see also |69|). The ob- 
tained eigenspectrum has a very rich structure. The quantized values of the conserved 
q— charges and the energy El depend on the integer rih and the real number defining 
the total ST(2,C) spin of the state, Eq. (3.16). In addition, they also depend on the 
"hidden" set of integers I = (£1, £2, ^2(l-2))- As a function of z/^, the charges form 
a family of trajectories in the moduli space q = (^2, <?3, <?l) labelled by integers rih 
and t. Each trajectory in the q— space induces a corresponding trajectory for the energy 
E L = E L {y h ]n hl t). The origin of these trajectories and the physical interpretation of 
the integers £ can be understood by solving the Schrodinger equation (3.7) within the 
semiclassical approach described in the next subsection. 

3.4 Semiclassical limit 

In the semiclassical approach |70|, we assume that the SX(2;C) spins h and h are large 
and apply the WKB methods to construct the asymptotic solution to (3.7). One might 
expect a priori that this approach could be applicable only for highly-excited states. 
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Nevertheless, as was demonstrated in |70|, the semi-classical formulae work with good 
accuracy throughout the whole spectrum. 

From the viewpoint of classical dynamics, the multi-gluon states (3.7) are describe 
by a chain of interacting particles 'living' on the two-dimensional b— plane |47,27|. The 
classical model inherits the complete integrability of the quantum noncompact spin mag- 
net. Its Hamiltonian and the integrals of motion are obtained from (3.9), (3.19) and 
(3.20) by replacing the momentum operators by the corresponding classical functions. 
Since the Hamiltonian (3.9) is given by the sum of holomorphic and antiholomorphic 
functions, from point of view of classical dynamics the model describes two copies of 
one-dimensional systems defined on the complex z— and z— lines. The solutions to the 
classical equations of motion have a rich structure and turn out to be intrinsically re- 
lated to the finite-gap solutions to the nonlinear equations |71,72|; namely, the classical 
trajectories have the form of plane waves propagating in the chain of L particles. Their 
explicit form in terms of the Riemann (9— functions was established in |47| by the methods 
of finite-gap theory |71,72|. 

In the semiclassical approach, the eigenfunctions in (3.9) have the standard WKB 
form, ^wkb^i, • • • , Zl) ~ exp(iS /h) where the Planck constant h = Iffel -1 ^ 2 is related 
to the lowest charge (3.21) and the action function So satisfies the Hamilton- Jacobi 
equations in the classical SL(2; C) spin chain. It turns out that the solutions to these 
equations are determined by the same spectral curve (2.27) as for the SX(2;R) spin 
chain. The charges q define the moduli of this curve and take arbitrary complex values 
in the classical model. Going over to the quantum model, we find that charges q are 
quantized. 

The quantization conditions for the charges q follow from the requirement that 
^wkb(2i, • • • , zl) has to be a single-valued function of As was shown in Refs. |27,70|, 
these conditions can be expressed in terms of the periods of the "action" differential over 
the canonical set of the a— and /3— cycles on the Riemann surface corresponding to the 
complex curve (2.27) 

Re <x> dxp(x) = tt £ 2 k-i i Re (p dxp(x) = tc £ 2k , (3.22) 

with k = 1, . . . , L — 2 and I = (£±, . . . , £21,-4) being the set of integers. The relations 
(3.22) define the system of 2(L— 2) real equations for (L—2) complex charges q^, ...,qL (we 
recall that the eigenvalues of the "lowest" charge q 2 are given by (3.21)). Their solution 
leads to the semiclassical expression for the eigenvalues of the conserved charges. In 
turn, the energy of the L— gluon states E L)q can be expressed as a function of q 3 , ...,qL- 
In the semiclassical approach, the corresponding expression is 

E { * s) = 4 In 2 (3.23) 



L 



+ 2ReJ^ 



fc=0 



V>(1 + iRe5 k + I Im6 k \) +ip(iReS k + | Im5 fe |) - 2^(1) 



where 5 k are roots of the polynomial t^u) defined in (2.24). 

The expression in Eq. (3.23) is similar to the energy of the SL(2,M) magnet in 
Eq. (2.28) although the properties of the two models are different. As was demonstrated 
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in |70|, the resulting semiclassical expressions for g 3 , and Ei are in good agreement 
with exact results |63,64|. A novel feature of the quantization conditions (3.22) is that 
they involve both the a— and /3— periods on the Riemann surface. This should be 
compared with the situation in the Heisenberg SL(2, M) magnet discussed in Section 2.6. 
There, the WKB quantization conditions involve only the a— cycles, Eq. (2.29), since 
the (3— cycles correspond to classically forbidden zones. For the 5X(2, C) magnet, the 
classical trajectories wrap over an arbitrary closed contour on the spectral curve (2.27) 
leading to (3.22). This fact allows us to explore the full modular group |73| of the complex 
curve (2.27) and explain the fine structure of the exact eigenspectrum of the SX(2;C) 
magnet. More details can be found in Ref. |70|. 

4 Concluding remarks 

In this review, we have described integrability symmetry in application to the deeply 
inelastic scattering in QCD. Due to space limitations, we did not review various impor- 
tant topics and we refer the interested reader to Ref. |53| for a comprehensive review on 
the subject. We would like to emphasize that integrability is not of a mere academic 
interest in QCD as it offers a powerful technique for solving important phenomenological 
problems such as finding the scale dependence of hadronic structure functions of higher 
twist and describing their high-energy (Regge) asymptotic behaviour. On theory side, 
the very fact that QCD evolution equations exhibit integrability property provides yet 
another indication that QCD possesses some hidden (integrable) structures waiting to 
be uncovered. 
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